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Abstract. In this paper I investigate what factors - both observational 
and physical - can change the measured slope of the observed 21cm HI 
power spectrum. The following effects can make the observed turbulence 
appear two dimensional rather than three dimensional: 1) if the turbu¬ 
lence is contained in a thin filament or slab; 2) if the medium has a high 
optical depth; and 3) if any method of observation or analysis is used 
which effectively limits the emission from the medium under study to a 
thin slab, for example, by analyzing an individual channel map. Straight¬ 
forward analysis of data can give misleading or incomplete results if these 
effects are not taken into account. 


1. Introduction 

The 21cm HI line has a column density power spectrum whose slopes are con¬ 
sistent with 2-Dimensional turbulence (a ~ —8/3) on large spatial scales (> 
0.01 pc) and narrow velocity ranges (Green 1993; Dickey & Crovisier 1983; 
Lazarian & Stanimirovic 2001; Dickey et al. 2001). The slope of this power 
spectrum is closer to a 3-D, Kolmogorov-like spectrum (a ~ —11/3) for wider 
velocity ranges but can still be significantly different than the Kolmogorov value. 
However, the electron density spectrum is consistent with Kolmogorov-like tur¬ 
bulence. This suggests a) that electrons and neutrals have different turbulent 
characteristics, or b) that there are effects which change the measured power 
spectrum slopes. 

Lazarian & Pogosyan (2000) have discussed how turbulent velocity may 
effect the observed power spectra of HI. However, there are other factors which 
Lazarian & Pogosyan did not discuss which may effect the power spectra. Among 
these are opacity and filamentary structures in the HI. 


2. Theory 

The power spectrum of the observed 21cm HI emission can be written as 

ps n {k) = c 2 n (k 2 + kt)~ a/2 ( 1 ) 

where k = , k Q = is the largest size scale of the turbulence and C 2 indicates 

the strength of the turbulence. Assuming that the observed HI intensity is 
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proportional to the column density it can be shown that 

(l{x)I(x + 5x)^ = C 2 {^N(x)N(x + 5x)^ (2) 


(l{x)I{x + 5x)^ 



( 3 ) 


where z and z' are the two lines of sight being compared. When the observations 
dictate that z Q «l Q then it can be shown that 


PS N (k X , ky) OC (kl + k 9 y+ k%) (i “ V2 


( 4 ) 


Thus, under common observing circumstances the column density power spec¬ 
trum can have an index that is one less than the index of the density power 
spectrum: ajv = a n — 1 for z a € (—oo, oo) and l < $(Z 0 ). 


3. Models 

Models of a turbulent cloud were developed by creating a Fourier Transformed 
data cube of the density. This was done by taking the real and imaginary parts at 
a given k from a Gaussian distribution with zero mean and a standard deviation 
given by the square root of the density power spectrum at that k. This data cube 
was then Fourier Transformed into a real density data cube. A simple radiative 
transfer scheme was then used to convert the data cube into the “observed” two 
dimensional HI brightness temperature. 

3.1. Limiting the Emission Depth 

To show that the above theory is correct, slices of different thicknesses were 
integrated from the same model to form the observed HI. The observed HI 
power spectra slopes as a function of thickness is shown in Figure 1. 

3.2. Opacity 

To study the effects of opacity, an average opacity from the front to the back of 
the data cube was used to define the detailed radiative transfer. As the opacity 
increases we can see from Figure 2 that the observed HI power spectra change 
slopes from rsj — 11/3 to rsj —3. 

3.3. Filamentary Structures 

The HI could also be confined to thin sheets or filaments. This view is supported 
by the images of 21cm HI emission that are being produced by the CGPS and 
SGPS surveys (see McClure-Griffiths 2002; and Knee 2002). Having the turbu¬ 
lence confined into filaments where vortices are stretched along the filament can 
also change the observed power spectral index. This can be realized in the model 
by changing k x to r/fc x where 0 < r] < 1 (i.e. making the turbulence anisotropic). 
In Figure 3 we compare the power spectra from a an isotropic model (?) = 1) 
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Figure 1. The power spectra of the same model when only the first 
30 pixels in z are integrated over and when the full data cube are 
integrated over. The model was created with a slope of —11/3. It is 
easily seen that any method of observing or a propagation effect that 
keeps the observer from seeing the entire turbulent cloud will change the 
observed power spectrum’s slope. The x-axis is the spacial frequency 
in 512 • pixels -1 and the y-axis is in arbitrary units. 



Opacity 


Figure 2. For the same model density cube this plot shows how the 
simulated observed HI power spectral index varies as the opacity (r) 
from the front to the back of the data cube is varied. This plot suggests 
that we may not be able to invert the observed power spectra back to 
the true power spectra since the function is not single-valued. 
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Figure 3. The structure functions of models of isotropic turbulence 
( r / = 1) and turbulence in filamentary structures (// = 0.1). Adding a 
value of two to the slope of the structure function and then multiplying 
by -1 gives the spectral index of the corresponding power spectrum. It 
is easily seen that for the filamentary case the power spectrum slope is 
—8/3 on the largest size scales. 

to a filamentary model ( r] = 0.1). It is easily seen that filamentary structure 
can produce a power spectrum with a slope of —8/3 on the largest size scales, 
1 > 1 filament where 1 filament is the thickness of the filament. 
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